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Abstract 

We investigate cooperative wireless relay networks in which the nodes can help each 
other in data transmission. We study different coding strategies in the single-source single- 
destination network with many relay nodes. Given the myriad of ways in which nodes can 
cooperate, there is a natural routing problem, i.e., determining an ordered set of nodes to 
relay the data from the source to the destination. We find that for a given route, the decode- 
and-forward strategy, which is an information theoretic cooperative coding strategy, achieves 
rates significantly higher than that achievable by the usual multi-hop coding strategy, which 
is a point-to-point non-cooperative coding strategy. We construct an algorithm to find an 
optimal route (in terms of rate maximizing) for the decode-and-forward strategy. Since the 
algorithm runs in factorial time in the worst case, we propose a heuristic algorithm that 
runs in polynomial time. The heuristic algorithm outputs an optimal route when the nodes 
transmit independent codewords. We implement these coding strategies using practical low 
density parity check codes to compare the performance of the strategies on different routes. 

1 Introduction 

Research in high data rate wireless systems has enabled applications to go wireless and become 
more interesting, e.g., wireless Internet access, mobile video conferencing and mobile TV on buses 
and trains. These applications would have been impossible without high rate wireless transmission 
links. As many wireless devices arc battery operated, power constraint is often imposed on them 
to make sure that they maintain a certain desired lifespan. In this paper, we investigate optimal 
routing problem to maximize the transmission rate in the wireless network where there is a power 
constraint on each node. 

The wireless channel is inherently broadcast, in that messages sent out by a node are heard 
by all nodes listening in the same frequency band and in communication range. This opens up 
opportunities for richer forms of cooperation among the wireless users/nodes. For example, rather 
than using point-to-point multi-hop routing (a direct adaptation from wired networks), where a 
node only transmits to the next node in the "route" , cooperative strategies, such as information 
theoretic relaying [3] and opportunistic routing ]4| , could be used. These richer forms of 
cooperation can lead to efficient distributed algorithms and can increase the end-to-end data 
rates. The gain from cooperation has been shown in information theoretic analyses [311] and 
demonstrated in practical implementations 015115]. 

We now briefly describe what we mean by these richer forms of cooperation, often using the 
term "coding" to highlight that our approach stems from information theory [TU]. Figs. 0a)-(c) 
depict wireless networks in which node 1 is the source, nodes 2 and 3 are relays, and node 4 
is the destination. In Fig. 0a), since every node can hear what node 1 transmits, the simplest 
strategy is for node 4 to directly decode from node 1, which we call the single-hop coding strategy 
(SH). However, when nodes 1 and 4 are situated far apart, signals from node 1 go through severe 
attenuation before they reach node 4. This is when relay nodes 2 and 3 can help. Referring to 
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Figure 1: Different coding strategies for multiple relay channels. 



Fig. [TJb), node 1 transmits to node 2. Node 2 fully decodes the data and re-transmits to node 
3. Node 3 does the same and relays the data to node 4. This is the well known multi-hop coding 
strategy (MH) . Although we can view relay nodes helping the source to transmit data as a form of 
cooperation, as far as decoding is concerned, SH and MH are still point-to-point strategies (a node 
only decodes from one node) and we categorize them as non-cooperative coding strategies. Taking 
a closer look at MH, we see that node 3 can hear and decode node l's transmission (although 
it is intended for node 2). This suggests a cooperative coding strategy, depicted in Fig. []Jc), in 
which node 3 decodes transmissions from nodes 1 and 2, and node 4 decodes transmissions from 
nodes 1-3. This cooperative way of encoding and decoding stems from an information theoretic 
approach and is termed the decode- and- forward coding strategy (DF) [UEIEI]. 

Regardless of whether MH or DF is used for data transmission, there is a sequence of nodes 
through which data flows. Kurose and Ross [IT] define a route as "the path taken by a datagram 
between source and destination" . The datagram "hops" from one node to the next node, capturing 
the scenario in which a node receives data only from a node behind and forwards data only to the 
node in front. However, in the cooperative coding paradigm, data does not flow from one node to 
another; rather it is from many to many with complex ways of cooperating. To describe the flow 
of information in these new modes of cooperation, we define a route as follows. 

Definition 1 The route taken by a packet from the source to the destination is an ordered set of 
nodes involved in encoding/transmitting and receiving/decoding of the packet. The sequence of the 
nodes in the route is determined by the order in which nodes ' transmit signals first depend on the 
packet. 

Remark 1 If a group of nodes transmits simultaneously, then they can be ordered arbitrarily 
within the group. For example, consider a four node network, in which node 1 first broadcasts the 
message, and then nodes 2 and 3 listen and simultaneously transmit to node 4- The route here 
can be described by {1, 2, 3, 4} or {1, 3, 2, 4}. 

Remark 2 Fig. [7] describes three coding strategies for the same four-node network. The route 
for SH in Fig. [TJa is {1,4}. The routes for MH and DF in Figs. {Tf) & \ljp respectively are both 
{1,2,3,4}. 

Given the myriad of ways in which nodes can cooperate, there is a natural routing problem 
in the cooperative coding paradigm. Furthermore, route selection directly affects the end-to-end 
data transmission rate. For DF, the current routing solutions for MH cannot be applied trivially. 
In this paper, we construct an algorithm to find an optimal route (in terms of maximizing rates) 
for DF. 

Our contributions in this paper are as follows. 

1. We show how much gain one can expect using DF, a cooperative coding strategy, over MH, 
a non-cooperative coding strategy, on the same route. 
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2. We construct an algorithm that finds rate maximizing routes for DF. 

3. We construct a heuristic algorithm that runs in polynomial time. We show that the heuristic 
algorithm finds an optimal route for DF when the nodes send independent codewords. 

4. We implement DF using low-density parity-check (LDPC) codes [12] [13]. Also, we show the 
performance of codes using different coding strategies and on different routes. 

This paper investigates cooperative coding and routing in the wireless network based on an 
information theoretic approach. A few idealized assumptions are made (e.g., infinite block length, 
unbounded communication range). Some of these assumptions are, however, relaxed in the simu- 
lations in Section [S] 

1.1 Related Work 

Communications in wireless networks has been progressing from MH to that using cooperative 
strategies. More research is being directed toward designing codes that are based on information 
theoretic cooperative coding strategies to harvest the gain in transmission rates predicted by 
information theory. Examples of codes based on cooperative coding strategies include DF-based 
Turbo codes [H] [IB] and LDPC codes SB E7J EH ]W\ for the single relay channel. It has been 
mentioned that some of these codes can be extended to the multiple relay channel 2\ |3J |5j [6 . 

In the past, link optimization (i.e., maximizing the transmission rate between node pairs) 
and route optimization were done separately. Routing was optimized after the links between the 
nodes had been established. Algorithms such as Bellman- Ford [2U[ Section 24.1] [21] and Dijkstra's 
algorithm [55J that assign costs to all links were used to find a route with the lowest cost from 
source to the destination. These ways of separating routing and coding are not optimal for MH 
or DF as the rates of the links change depending on which route is chosen. Realizing the inter- 
dependency between links and routes, it has been suggested that links and routes be jointly 
optimized [531 HH HS1 US] ■ This gives rise to cross-layering [57] in the OSI model. However, in 
these joint routing and coding work, data transmission from the source to destination is still based 
on MH. Routing algorithms that are optimized for MH might not be suitable for DF. 

In Ad hoc On-demand Distance Vector Routing (AODV) [55] and Dynamic Source Routing 
(DSR) [29], the source node broadcasts a route discovery packet. Neighboring nodes receive and 
re-broadcast the packet. When the destination receives the packet, a route is formed by tracing 
the path that the packet took. These routing algorithms minimize the transmission delay but 
might not optimize the transmission rate. 

In Extremely Opportunistic Routing (ExOR) [I] , a node broadcasts its data to a set of potential 
relays. Nodes in this set transmit acknowledgments and then selected nodes forward the data. 
Though ExOR does not have predefined routes, MH is used on the effective route taken by a 
packet. 

As far as we know, routing algorithms for cooperative coding have not been investigated. In 
this paper, we propose algorithms to find optimal routes for DF-based codes in the multiple relay 
channel. Our work complements code design by finding the best route (rate maximizing) on 
which the codes can be used. As previous work focused on cooperative coding for the single relay 
channel, in this paper, we implement DF-based LDPC codes on the multiple relay channel. We 
then compare the transmission rate of different coding strategies on different routes. 

We focus on the multiple relay channel [5J [5] [S] [S] , which is a single-source single-destination 
network, as a first step towards understanding general multiple-source multiple-destination net- 
works. We study DF because it is one of the "more implementable" information theoretic coding 
strategies Ql HI EH El H3 US] • 



3 



2 Motivating Cooperation 
2.1 Network Model 

We consider a D-node network S = {1,2,3..., D — 1,D} with one source (node 1) and one 
destination (node D). Node «,V« £ S, cither transmits at fixed average power Pj or turns off. 
We use the standard path loss model for signal propagation. The received power at node t from 
node i is given by Pu — nd^Pi, where da is the distance between nodes i and t, r\ is the path 
loss exponent {r\ > 2 with equality for free space transmission), and k is a positive constant. The 
receiver at node t is subject to thermal ambient noise of power N t . We assume duplex nodes, i.e., 
nodes can transmit and receive simultaneously. We assume that all nodes have the same noise 
variance. 

Given this network model, we investigate how nodes can cooperatively send messages from the 
source to destination. Wc study and compare several coding strategies. 

Remark 3 We consider single-flow networks. This is the first step in understanding a more 
complicated problem of multiple flows. The relevance of our work in multiple-flow networks is as 
follows: 

1. In a multiple-flow networks where each flow uses an allocated orthogonal channel, the rate 
of each flow can be optimized in respective channel using the algorithm derived in this paper. 

2. In a multiple-flow network with existing flows, if we wish to add a new flow, the algorithm 
in this paper finds an optimal route for the new flow. Note that adding a new flow might 
affect existing flows. We can restrict the transmit power of nodes in the new flow to control 
the interference introduced. 



In SH, the source directly transmits data to the destination. The signal-to-noise ratio (SNR) at 
the destination, node D, is given by 7sh(-D) = PujNp 1 . The Shannon capacity of this SH link is 
i?SH = 5 log (1 + 7sh(-D)) . This rate depends on the source-destination distance and can be poor 
if the source and destination are situated far away from each other (because of signal attenuation) . 

Remark 4 We assume that nodes that do not participate in relaying the data for a source- 
destination pair do not cause interference. Another way to account for the external noise is to 
include it in the receiver noise. 



In MH, we make use of the relays to aid the transmission from the source to the destination. The 
source simply transmits to the next relay. The first relay decodes the message and re-transmits 
it to the second relay, and so on until the destination. This can improve the transmission rate 
if the attenuation from the source/relay to the next relay is reduced as compared to that in SH. 
However, since all relays transmit simultaneously, there exists interference, beside noise, at the 
receiver. 

In the rest of this paper, we denote a route by M = {mi, m 2 , . . . , m\ M \}. We define the 
set of all possible routes from the source (node 1) to the destination (node D) by H(S) = 

< {mi, mi, . . . , m\ M \\ : m 2 , . . . , m\M\-i are an possible selections and permutations of the relays 
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Single-Hop Coding Strategy (SH) 



2.3 Multi-Hop Coding Strategy (MH) 
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Since all relays and the destination must fully decode the messages, the transmission rate from 
the source to the destination using route A4 is 

Rmk{M) = min -log(l+j M n(m,M)). (2) 

mEA"f\{rai} Z 

We term Rmh(-M) the rate supported by the route M. using MH. The maximum rate using by 
MH, optimized over all possible routes, is 

RZn= max Rmk(M). (3) 
Men(s) 

We note that there may exist more than one route that support this maximum rate. 
2.4 Decode-and-Forward Cooperative Coding (DF) 

Using DF [U[2j[3], each decoder decodes transmissions from all nodes behind. E.g., the third node 
in the route decodes the transmissions from the first and the second node. So, the third node 
decodes each source message using two blocks of received codewords. In addition, assuming that 
the nodes in front decode the messages correctly, a node knows what they transmit and hence it 
cancels the interference from these nodes. It has been shown that in order to maximize the DF 
rate on route M, node m, transmits X m% = Y/j^j+i y / a mi m j P mi U m:j , for < Yjj^i+i «m, mj < 1, 
= 1, . . . , \M.\ — 1. U m . are independent Gaussian random variables with unit variance, {a^ \j = 
i + 1, . . . , | M. |} are the power splits of node i, allocating portions of its transmit power to transmit 
independent sub-codewords Uj. Doing this the SNR of node m t in route M is 




7df(to 4 , 

*0 = ^ X X \ a m ^P mtmt . (4) 



Using the route Ai, DF can achieve rates up to 



Rdf(M) = max min - log (1 + 7dp(to, M)) , (5) 

°ij meJM\{mi) 2 

and the maximum rate using by DF is 

R™*= max Rbf(M). (6) 

Men(s) 

Definition 2 We define the reception rate of node m in route as R m {M) — h log (1 + 7£>i?(m, M.)) . 
It is the rate at which node m can fully decode the messages. The same concept applies to MH. 

Remark 5 In practice, a relay in the route might decode a message wrongly, and hence forward 
the wrong message. When this happens, the nodes behind, when trying to cancel the co-channel in- 
terference introduced by this relay, will introduce more noise at their decoders. While this scenario 
is not captured in (j4|) ; we allow imperfect interference cancellation in our simulations (Section^. 



2.5 Comparing the Strategies 

It is easy to see that, for any chosen route M. with four nodes or more, 7df("t-, M) > 7mh(w, M), Vto £ 
A4. Also, we can show that for any M. £ 11(5), 

Ruf(M) = Rmh(M) = R S h, far|A<| = 2, (7a) 

Rbf{M) > Rmh(M), for \M\ = 3, (7b) 

R DF (M) > R U h(M), for \M\ > 4, (7c) 
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Figure 2: Ratio of average transmission rate SH, MH, and DF versus \M\ for two cases: (1) 
d w = 10 and (2) d 1D = \M\ - 1. 

and i?g| x > -R^iff — ^SH- 

However, it is not clear how much again, on average, we can expect using DF compared to 
MH and SH. Now, we compare the rates of SH, MH, and DF for randomly generated routes of 
different lengths in a line topology. We consider two cases: (1) dio = 10 and (2) dm = \M- \ — 1- 
Then |A^| — 2 nodes are randomly placed along the straight line joining nodes 1 and D. Note that 
in case 1, node density increases with the number of nodes while in case 2, the average adjacent 
node spacing is constant for all \M\. We set Pj = JVj = 1 for all transmitters and receivers, and 
K = 1, rj = 2. For each randomly generated route, we calculate the transmission rate using SH, MH, 
and DF. Here we restrict the nodes to transmit independent codewords for easier optimization, 
i.e., we set onj = 1, Vi, Vj = i + 1 and Oy = 0, Vj ^ i + 1. 

Fig [5] shows the results for cases 1 and 2. For a route of 25 nodes, the DF rate is roughly two 
orders of magnitude higher than that of SH and MH for both cases. Moreover, as we increase the 
number of nodes in the route, the gain of DF over MH/SH increases for both cases. 

We note that if the nodes are allowed to send arbitrarily correlated codewords, the DF rate 
can be higher. For example, consider the route M = {(0, 0), (0.5, 0), (2, 0), (3, 0), (4, 0)}. With 
independent codewords, 

Rdf(-M) I Rmh(-M) = 2.95, but with arbitrarily correlated codewords, 
R BF (M)/R M n(M) = 4.40. 

3 The Optimal Routing Problem 

We define the optimal route set for DF as 

Qdf = {ME U(S) : R DF (M) - R^} (8) 

where H(S) is a set of all possible routes from the source to the destination. We define the optimal 
route set because the rate maximizing route may not be unique. Then the optimal DF routing 
problem is 

Find at least one -M^p & Qdf and R BF «) . 
The optimal route set and routing problem for MH are similarly defined. 
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Finding optimal routes for MH and DF by brute force is hard, as it involves testing all routes 
in II(£>). In Sections |4] and El we construct an algorithm that finds potentially without 

having to test all routes in 11(5). However, this algorithm runs in factorial time in the worst case. 
Hence, in Section [71 we proposed a heuristic algorithm, which runs in polynomial time. 

4 The Nearest Neighbor Algorithm 

Now, we present an algorithm to find an optimal route for DF. In the section, we assume that 
nodes use independent codewords, i.e., we set ay = l,Vi,V) = i + 1 and ay = 0, Vj ' =^ i + 1. 

Remark 6 Although the theorems in this section are proven assuming that the nodes send in- 
dependent codewords, we can also show that they hold even when the nodes transmit arbitrarily 
correlated codewords j30f . 

Remark 7 We consider independent codewords in this section because coherent combining is 
practically infeasible. When the nodes are operating in the GHz range, it is difficult, if not im- 
possible, to synchronize the carriers to nanosecond accuracy. Furthermore, even if we have very 
precise clocks, coherent combining is still unlikely in a multiple-node network. For example, in a 
four-node route, even if we manage to synchronize nodes 1 and 2 to allow coherent combining at 
node 3, we might not be able to ensure that they will coherently combine at node 4- 

First, we define the nearest neighbor with respect to a route. 
Definition 3 Node i ^ M. is a nearest neighbor with respect to the route A4 iff 

P mi > Pmj, Vm e M,yj eS \ (M U {i}). (9) 

Note that nearest neighbor might not be unique. Now, we describe the nearest neighbor algorithm 
(NNA). 

Algorithm 1 (NNA) 

1. Initialize M. = {mi}, where mi = 1. 

2. If there exists a unique nearest neighbor i* with respect to the current route M, we append 
i* to the current route: Ai *— M. U {i*}. Else, the NNA terminates prematurely. Since M. 
is an ordered set, the notation A U B means appending ordered set B to the end of ordered 
set A. 

3. Step^is repeated until the destination node, node D, is added into M.. 

The algorithm is said to terminate normally if node D is added to the route. Otherwise, the 
algorithm is said to terminate prematurely. If the NNA terminates normally, we have the following 
theorem. 

Theorem 1 Consider a multiple node wireless network with one source and one destination. If 
the NNA terminates normally, then the NNA route is optimal for DF. 

To prove Theorem [TJ we need the following lemmas. 

Lemma 1 When we add the unique nearest neighbor, node a* , to route M., the rate supported by 
the new route M.i — M. U {mijfj = a*} is greater or equal than the rate supported by the route 
formed by adding any other node to M.. Mathematically, 

Rdf{M U {a*}) > R DF (M U {b}), V6 e S \ (M U {a*}). (10) 



7 



Proof: [Proof for Lemma [T] Considering M2, the reception rate of node 



m\ M2 



= b is 



R b (M U {b}) = - log 



\M\ 
i=l 



(11) 



and the reception rate of the node TOi^i = a* in route Mi is 



i? Q .(XU{a*}) 



■log 



\M\ 

i + n-}J2 p ™^ 



(12) 



Clearly, if P ma * > P m b, Vra G M with at least one inequality, R a *{Mi) > Rb{M.2)- Hence 
Rvf(Mi) > i?DF(-A^2)- Hence we have Lemma[TJ 

We have proven that at any point of time during route construction, in order to maximize the 
rate supported by the route, we must choose the nearest neighbor (assuming it exists). Next, we 
show that choosing the nearest neighbor will not harm the rate supported by the route even when 
more nodes are added. 

Lemma 2 Let A4 = {a^ai,, ■ ■ ■ , a *M\\ be a row ^ e formed by adding the nearest neighbor one by 
one starting from the source. Now, arbitrarily add K nodes to M. The first node b\ is not 
a nearest neighbor and the rest may or may not be nearest neighbors. In other words, M.\ = 



'2' ' ■ ■ >"|X| 



61, 62, ■ ■ ■ , bx}, where b\ is not a nearest neighbor to M.. We can always replace 



b\ by the nearest neighbor a? M , +1 (assuming it exists) to obtain 



{a*, . . . , o.* M \,a* M \ +1 ,bi, . . . , b K -i}, 

l f a \M\+i i {h,...,b K -i}, 
{a*, . . . , a*M\i a \M\+i>bi> ' ■ ' ' °k-i,bk+i, 
a *\M\+i = b k, for some b k G {61, . 



M 2 = 



■,ok}, 

,OK-l}, 



(13) 



where Rdf{M2) > Rdf(Mi 



Proof: [Proof for Lemma [2] For both cases in (fT3"|) . the reception rates for the first \M\ nodes 
in both M.\ and M2 remain the same as each of them decodes from the same nodes behind the 
route. In equations, 

Ra t {Mi) = Rar {Mi), W = 2,3,...,|A4|. (14) 

From LemmaCQ R a * M X {M%) > Rb 1 {Mi). 

Now, we study the case when o,* M ^ +1 ^ {b\, . . . , bx-i}- For nodes {b\, . . . , &rt-i} in M2, with 
an additional node behind, i.e., o* M ^ +1 , the reception rates of these nodes are higher than the 
same nodes in Mi: 

R bi {M 2 ) >R bi {Mi), Vi=l,2,...,K-l. (15) 

Now, we study the case when a* M ^ +1 = b^ for some b^ G {&i, . . . , bK-i}- Similar to the first 
case, with an additional transmitting node a* M ^ +1 , the nodes {&i, . . . , in M2 have higher 

reception rates compared to those in Mi, i.e., 



R bi {M 2 ) >R bi {Mi), V* = l,2,. 



1. 



(16) 



For nodes {b^+i, ■ ■ ■ , 6 if}, there is no change in the reception rate because each of them has exactly 
the same nodes behind them in both M2 and M\. So, 



Rb i (M 2 ) =R bi {Mi), Vi = k + l,k + 2,...,K. 



(17) 



Lemma[2]follows by noting that Rf>f(M 2 ) = m ^ n ieM 2 \{a*} Ri{M2) > niin ieMl y a .) Ri(Mi) = 
Rdf(M 1 ). 
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Lemma 3 For a route that contains all nearest neighbors, the supported rate is always higher or 
equal to any route, of the same length, with one or more non-nearest neighbors in it. 

Proof: [Proof for Lemma [3J Lemma [3J can be proven by applying Lemma [3] recursively until 
the entire set is replaced by nearest neighbor nodes. 

Now we consider routes of different lengths but that end on the same node. 

Lemma 4 Consider a route Mi, where node 1 is the source, and node m\M x \ — D is the desti- 
nation, that is 

Mi = {m* 1 ,m 2 ,...,m\ Ml \}. (18) 

Here, one or more nodes in {m 2 , ■ ■ ■ , TO|a/i|} are not nearest neighbors. The following route where 
all nodes are added according to the NNA (assuming that it does not terminate prematurely) 
supports rate as good or higher than that supported by Mi- 

M 2 = {m\,ml,...,m* M ^}, (19) 

where m^ M ^ = D and \M\\ not necessarily equals \M 2 \. In other words, Rdf(M 2 ) > Rdf(Mi). 

Proof: [Proof for Lemma U First of all, we consider the case \M\\ — \M 2 \. The results follows 
immediately from Lemma [3J Second, we consider |Mi| > \M 2 \. We consider first \M 2 \ nodes in 
Mi, i.e., M[ = {m\,m 2 , . . .,m\M 2 \}- Then, 

Rdf(M 2 ) > RMM[) > Rdf(Mi). (20) 

The first inequality is obtained by applying Lemma [3J M 2 and M\ are of the same length. The 
former is formed using the NNA while the latter is not. The second inequality can be argued as 
follows. The first \M 2 \ nodes in both routes M{ and Mi are identical. Hence the reception rates 
are the same. However, there are additional nodes in Mi whose reception rate might be lower 
than Rdf(M' 1 ). Hence, the rate supported by M[ can only be higher than that of Mi. 

Lastly consider IM2I > \M\\. We replace the transmitting nodes in Mi with nearest neighbors 
and obtain 

M 3 = {ml,m 2 , . . . ,rri\ Ml \-n m \Mi\ = D }- ( 21 ) 
Note that roi^i might not be the nearest neighbor. Clearly, using Lemma [H Rdf{Ms) > 
Rf>f{Mi). Now, 

Rdf(M 2 ) 

= min{i? m «(X 2 ), . . .,R m j Miii (M 2 ),R m j Mi] (M 2 ), 

. . . , R mU ^{M 2 ),R m ^ (M 2 ) = R D (M 2 )}, 

> mm{R rn ,(M 3 ), Rm*^ (M 3 ), R D (M 3 )}, 

= Rdf{M 3 ) > R OF (Mi). 

The inequality in (|22b[) is because in M 2 , { m J Ml \, ■ ■ ■ , m *M 2 \-i^ are a( ided to {m\, 
before D. A necessary condition for this is 

Pmn > PrnD, Vm £ {m*, . . . , m* M ^_ 1 }, 

Vn e {m[_ Mi| ,...,mj t _ M2| _ 1 }, (23) 

with at least one strict inequality for each n. Hence, R n {M 2 ) > Ro{M 3 ),^n £ {m^ M ^ , . . . , mf M \_y\. 
With additional nodes transmitting to D in M 2 , Rd{M 2 ) > Rd(M 3 ). Hence, we have Lemma|H 
Proof: [Proof for Theorem [T] From Lemma [H we know that if the NNA terminates normally, 
the route (from the source to the destination) formed using the NNA can support transmission 
rates as high as any other route. In other words, the NNA finds a route that supports the highest 
rate achievable by DF. Theorem [1] follows. 

Remark 8 We note the NNA terminates normally if and only if a unique nearest neighbor exists 
at each step. In the next section, we extend the NNA to an algorithm which terminates normally 
given any network topology. 



(22a) 
(22b) 
(22c) 

•••' m t«i|-i> 
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5 The Nearest Neighbor Set Algorithm 



In this section, we modify the NNA so that it terminates normally in any multiple node wireless 
network with a single source and a single destination. We term this algorithm the nearest neighbor 
set algorithm (NNSA). First, we define the nearest neighbor set. 

Definition 4 The nearest neighbor setAf — {n±, n%, ■ • ■ , fiijvi} with respect to route A4 = {mi,m2, 
. . . * s defined as the smallest set AT where each n G Af C S \ M. satisfies the following 

condition. 

Pmn>P ma , VmeM,VaeS\(MUAf), (24) 
with at least one strict inequality for every pair of (n, a) G {(n, a) : n G Af, a G S \ (A4 U A/")}. 

Now we describe the NNSA. 
Algorithm 2 (NNSA) 

1. Starting with the source node, we have A4 = {1}. 

2. Find the nearest neighbor set Af. The original route AA branches out to \Af\ new routes as 
follows: 

M t <- M U i=l,...,\Af\. (25) 

3. For each new route in (|25p . step\^is repeated until the destination is added to all routes. 

When the algorithm terminates, we end up with many routes from the source to the destination. 
We term these routes NNSA candidates. We calculate the supported rate of each candidate and 
choose the one which gives the highest supported rate. The following theorem says that any NNSA 
candidate that gives the highest supported rate is an optimal route for DF. 

Theorem 2 Consider a single-source single- destination multiple node wireless network. The 
NNSA candidate routes that give the highest supported rate are optimal for DF. 

Proof: [Sketch of proof for Theorem [2] Using the technique used in the proof of Theorem [TJ 
we can show that adding a node that does not belong to the nearest neighbor set can only 
be suboptimal. We can always replace that node with one from the nearest neighbor set and 
obtain an equal or higher rate. In other words, we can show that the supported rate of A4\ = 
{1, m.2, . . . , mijifj—ij D}, with one or more nodes in {m.2, . . . , toij^i} not from the nearest neighbor 
set, is lower or equal to the supported rate of M.^ — {1, m^, ■ • • , m *M 2 \-i> w here all nodes in 
M.-2 are added according to the NNSA. In other words, i?DF(A^2) > Rdf(-Mi). 

The NNSA finds all possible routes for which every node is added from the nearest neighbor 
set. Hence one or more of the NNSA candidates must achieve the highest rate achievable by DF. 
This gives us Theorem [2j 

Remark 9 We can show that a shortest optimal route, defined as some M^p R G Qdf, s.t. 
\-M.df R \ ^ |A4|, VA4 G Qdf, is contained in one of the NNSA candidates that supports R^p x . 

Remark 10 The NNSA might output optimal routes that include more nodes from the network 
unnecessarily. In other words, shorter optimal routes exist. However, from Remark\Q we can find 
the shortest optimal route by pruning the optimal routes output by the NNSA . 
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6 Complexity of NNSA 



With the NNSA, we can now search for the optimal route in the NNSA candidate set, as compared 
to searching in H(S) using brute force. The number of candidates determines the number of routes 
whose rate we need to calculate to find optimal routes for DF. We note that the size of the NNSA 
candidate set might still, in the worst case, equal |II(<S)|. Using brute force, the number of 

permutations we need to check is |II(<S)| = 1 + ( D ~ 2 ) + ( D ~ 2 ) H (d-I) = °(( D ~ 1 )0> where 

D is the total number of nodes in the network and (^) = ( n Jtyx™n-k-'\) x » x 1 • 

We ran the NNSA on 10000 randomly generated networks with a varying number of nodes 
uniformly distributed in a lmxlm square area. The source, relays and destination were randomly 
assigned. On average, half of the NNSA candidate set sizes were less than 0.715% of the size of 
II(6>) for the 8-node channel and less than 0.253% for the 11-node channel. 

We note that the average size of the NNSA candidate set does grow factorially with the number 
of nodes. However this does increase the range of finite size networks for which we can find optimal 
routes. Furthermore, the NNSA provides insights for designing heuristic algorithms to find good 
routes for DF based codes. In the next section, we propose heuristic algorithms which find routes 
in polynomial time. 

Remark 11 The NNSA builds routes from the source. Relays are added to the route regardless 
of where the destination is. We will use this observation in designing heuristic algorithms. 



7 Heuristic Algorithms 

In the NNSA, the optimal route is constructed by adding the "next hop" node one by one to the 
partial route. The node to be added is from the nearest neighbor set. If the nearest neighbor set 
contains more than one node, the current route branches to more than one route, leading to a 
possibly large NNSA candidate set size. 

To avoid this, we consider a heuristic approach that starts from the source node and repeatedly 
adds only one "good" candidate from the nearest neighbor set until the destination is reached. For 
the choice of the next hop node, we consider the node which receives the largest sum of received 
power from all the node in the existing route. We call this the maximum sum-of-received-power 
algorithm (MSPA). By choosing only one node to be added to the partial route, we prevent the 
algorithm from branching out to multiple routes. This heuristic approach yields only one route, 
regardless of the network size. We now explicitly describe the MSPA. 

Algorithm 3 (MSPA) 

1. Start with the source node: M. = {1}. 

2. For every node t £ S\A4, find the sum of received power from all nodes in M. to t, X]ieA<f 

3. Let a* be any node with the highest sum of received power, i.e., J2ieM ^ ia * — J2jeM Pjt^t S 
S\M. Append node a* to the route: M <— M U {a*}. 

4- Repeat steps until the destination is added to the route. 

Remark 12 Assuming that the value of the previous sum-of-received-power computations are 
cached, the complexity of step [H in MSPA is O(D) because there are at most (D — 1) nodes 
not in the route. The complexity of the comparisons in step\^is 0{D). Steps HH2| are repeated at 
most (D — 1) times, giving a worst case complexity of the MSPA ofO(D 2 ). Recall that D = \S\. 

It turns out that the MSPA is optimal if the nodes are restricted to sending independent 
codewords, as proven in the following theorem. 

Theorem 3 In a single-source single- destination multiple node wireless network in which the 
nodes send independent codewords, the MSPA route is optimal for DF. 
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Proof: [Proof for Theorem [3] Consider an optimal route 
Mi = {ml, m%, ■ ■ ■ , m* k , m^ +1 , . . . , m* M ^}. Suppose that the first k nodes of the MSPA route 
are the same as this optimal route but the (k + l)-th node is different, i.e., the MSPA route is 
M2 = {m\, ml, . . . , ml, a, . . . } where a ^ ml +1 . 

Since a is added to the route by MSPA, a necessary condition is 

Ra(M 2 )>R m i +1 (Mi). (26) 

Now, consider the case where a ^ m*,Mi = k + 2, . . . ,\M\. We add a to Mi and obtain 
M 3 = {ml, ml, . . . ,m* k ,a,m* k+1 , . . . ,m* M |}. Then, 

Rm*{M 3 ) =R mt (M 1 ), i = 2,...,k (27a) 
R a (M 3 ) > Rm l+l (Mi) (27b) 
Rm t (M 3 ) >R mJ (M 1 ), i = k + l,...,\M\. (27c) 

So, R DF (M 3 ) > Rbf(Mx). 

Suppose a — m* , for some n G {k + 2, . . . , |jM|}. We swap the position of a and obtain 
Ma = {ml,m* 2 , . . . ,m* k ,a,m* k+1 , . . . , m*^, m* +x , rn]* M |}. Then, 

Rml (Mi) =Rm*{Mi), i = 2,...,k,n+l,...,\M\ (28a) 
RaiMi) > R m i +1 (M!) (28b) 
Rm*{Mi) >R m *(Mi), i = k + l,...,n-l. (28c) 

So, R DF (M±) > Rbf(Mx). 

In summary, we choose an optimal route. Starting from the second node, we compare the 
optimal route with the MSPA route. If the nodes are different, we insert (or swap, if the node is in 
the optimal route but at a different position) the node in the MSPA route to the optimal route to 
obtain a new optimal route. Repeat this by comparing the new optimal route to the MSPA route 
and changing the first node that differ until the MSPA is contained in an optimal route. Then, 
we have Theorem [31 

Remark 13 However, unlike the NNA and the NNSA, the MSPA does not output an optimal 
route when the nodes are allowed to send arbitrarily correlated codewords. Consider a four-node 
network with node coordinates 1(0,0), 2(04 18,0), 3(0.209,0.6755), and 4(0.995,0). Assume P % = 
l,Ni = 1,k = 1,7] = 2. The MSPA route is Mi = {1,2,4}. The NNSA outputs Mi and 
Mi = {1,2,3,4}. It is easy to compute that R DF (Mi) = 1.30826 and R DF {M 2 ) = 1.31576. 



8 DF with LDPC Codes 

In the previous section, we computed achievable rates of different strategies and routes based on 
an information theoretic approach. In this section, we compare the different strategies and routes 
in a line network using practical low density parity check (LDPC) codes [12] [13] with incremental 
redundancy. The aims of this section are: 

1. To illustrate that DF on the multiple relay channel is implementable. 

2. To demonstrate that DF performs better than MH under certain network topologies. 

3. To show that routing backward (away from the destination) can be good in DF. 

4. To demonstrate that the NNSA route performs better than other routes using DF. 
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Figure 3: Network topology. 
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Figure 4: Performance (information bit error rate (BER) versus transmit SNR) of different strate- 
gies on different routes in a 4-node network. 

From Fig. [5J we see that for a given route, DF performs better than MH. An interesting 
observation is that routing backward helps in DF but not MH. We find that the NNSA route 
(which is also the MSPA route), i.e., {1,2,3,4}, achieves the lowest BER compared to other 
routes using DF. 

Remark 14 We note that the total transmit energy differs depending on the length of the route. 
One might argue that route {1,4}, though having a higher BER, is better as only 1/3 power is 
consumed compared to route {1,2,3,4}. However, we stress that this paper finds a route that max- 
imizes the transmission rate, given that each node must transmit within a given power constraint. 
Whether or not the node transmits, it is not important in the route comparison. 

Remark 15 We plotted BER versus SNR in Fig. [7J If the maximum raw channel data rate 
(in bps) is $ m(lI , then the throughput is ^> = (1 — PER) ■ \t max > where PER is the packet error 
rate and depends on the BER and the packet size. In simulations, we found that packet error rate 
(PER) had the same behavior as BER. 



9 Concluding Remarks 

We first showed that DF gives a significant transmission rate gain over MH, for an arbitrary 
route in the wireless network. We presented an algorithm, the nearest neighbor set algorithm 
(NNSA), to find optimal routes, which maximize the rates achievable by DF. As this algorithm, 
in the worst case, runs in factorial time, we designed a heuristic algorithm, the maximum sum- 
of-received-power algorithm (MSPA), that runs in polynomial time. We showed that the MSPA 
finds an optimal route when the nodes can only send independent codewords. However, unlike 
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the NNA and the NNSA, the MSPA does not find an optimal route when the nodes are allowed 
to send arbitrarily correlated codewords. We implemented DF on practical networks using LDPC 
codes with incremental redundancy to compare different routes. 

We would like to highlight that for a given route, the choice of coding strategies, be it MH or 
DF, does not affect the spatial re-use of the system. In both strategies, the nodes in the route 
(except the destination) transmit at the same power level. The difference lies in how the nodes 
decode the data. 

We also note that there are some practical problems in implementing DF in large networks. 
First, since real-world systems have finite transmit power constraints, nodes have a finite commu- 
nication range, beyond which they cannot be heard. Second, large networks may be distributed 
over a wide area and cooperation over large distances may not be feasible. The solution is to 
partition the network into clusters and perform cooperative coding locally in the cluster, e.g., local 
DF, in which only nodes in a cluster cooperate with each other. 
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